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SUMMARY 


The 'objective  of  this’project'was  to  establish  computer  simulation  as  a  tool 
for  studying  scattering  from  randomly  rough  surfaces.  The  proposed 
methods  to  be  investigated  included  the  standard  moment  method  and  the 
integral  equation  method  (IEM). "  The  project  began  by  investigating  the 

u*>a  •  £  *  ’  ■*  J  -for  **■ 

minimum  size  e{  the  surface  patch  necessary  to  obtain  the  scattering  cross 
section  per  unit  area,  which  should  be  the  same  as  obtained  from  an  infinite 
surface  patch  when  the  method  used  was  the  standard  moment  method.  -It  <-■ 
^  _was  fgund-that  five  points  per  wavelength,  or  per  surface  correlation  length 
(whichever  was  the  smaller),  were  needed  and  the  minimum  surface  patch 
should  be  at  least  seven  electric  wavelengths,  or  surface  correlation  lengths 
(whichever  was  the  larger).  The  minimum  number  of  surface  patches  to  be 
averaged  to  obtain  a  satisfactory  estimate  of  the  scattering  coefficient  was 
25.  This  means  that  the  smallest  number  of  unknowns  to  be  evaluated  was 


for  a  scaler  problem  and  2  times  vector  problem.  To  determine 

this  many  unknowns  for  a  scalar  problem  one  needs  to  invert  a  matrix  size 
of  $5}  by  35^  in  the  moment  method .  In  general,  the  size  of  the  problem 
would  be  larger  and  hence  the  integral  equation  method  was  investigated. 

^^as-found-fchat  the  IEM  was  a  better  method  than  the  moment  method  in  that 
* 

it  required  much  less  memory.  It  required  a  much  greater  computational 
time,  heweve^  For  this  reason  another  method  called  the  frequency  domain 
method  (FDM)  was  considered  towards  the  end  of  the  project.  The  FDM 
was  the  preferred  method^as  compared  with  the  moment  method  and  the 
IEM^iirtwo-dimensional  surface  scattering  problems  because  it  was  more 
efficient  in  both  memory  and  computational  time.  Due  to  time  limitation. 
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extensions  of  the  FDM  to  three  dimensional  scattering  problems  are  yet  to  be 

.  I  .  /  i 

done.  ,  . :  / 

t 

f  , 

1.  INTRODUCTION 

This  project  was  started  in  August  1987  but  the  binding  was  discontinued  in 
February  1988.  It  was  restarted  in  January  1989  for  a  period  of  one  year. 
The  proposed  studies  were  completed  and,  in  addition,  a  new  approach,  the 
frequency  domain  method,  was  also  investigated.  The  chronological 
progress  of  the  project  is  given  in  the  next  paragraph  and  accomplished 
studies  are  given  in  the  next  section  with  reference  to  the  reports  delivered 
during  the  course  of  the  project.  The  work  done  over  the  last  quarter  is 
given  in  Section  3.  A  conclusion  is  given  in  Section  4. 

The  aim  of  this  project  was  to  establish  computer  simulation  of  radar 
scattering  as  a  tool  for  investigating  surface  scattering  phenomena  by 
defining  the  scope  of  the  problem  and  the  technique  to  be  used.  At  the  time 
the  proposal  was  being  prepared  the  standard  method  for  simulating  surface 
scattering  was  the  method  of  moment.  The  use  of  the  moment  method 
requires  the  inversion  of  a  matrix  of  size  proportional  to  the  square  of  the 
number  of  unknowns.  In  sun  ...  scattering  the  number  of  unknowns  is 
equal  to  the  number  of  surface  ^  its  needed  to  define  the  surface  patch 
from  which  we  want  to  compute  scattering.  Thus,  for  a  surface  with  100  by 
100  points  we  have  a  matrix  of  size  104  x  104,  which  is  clearly  too  large  for 
today's  supercomputer.  Hence,  an  immediate  interest  is  to  determine  how 
small  a  surface  patch  we  need  to  simulate  the  scattering  cross  section  per  unit 


area  from  an  infinite  patch.  The  immediate  goal  was  to  find  the  answers  to 
the  following: 

(a)  What  was  the  minimum  size  of  the  surface  patch  in  terms  of  the 
number  of  surface  correlation  lengths  or  the  electrical  wavelengths-- 
whichever  was  larger? 

(b)  What  was  the  minimum  number  of  points  per  electrical  wavelength 
or  surface  correlation  length-whichever  was  smaller? 

In  the  event  that  the  minimum  requirement  for  the  number  of  unknowns 
came  out  to  be  larger  than  what  was  suitable  for  the  moment  method,  it  was 
proposed  that  the  integral  equation  method  (IEM)  should  be  investigated.  It 
was  demonstrated  that  EEM  was  a  viable  method  which  could  be  handled  by 
to-day’s  computer.  However,  it  took  a  lot  of  computer  time.  Finally,  a 
frequency  domain  method  was  investigated  and  demonstrated  to  be  efficient 
in  both  memory  and  computational  time  for  a  two-dimensional  scattering 
problem.  The  theoretical  basis  for  its  extension  was  developed  during  the 
last  period  of  this  project  and  is  given  in  Section  3. 

2.  STUDIES  ACCOMPLISHED 

Surface  scattering  simulation  begins  with  a  computer-generated,  statistically 
known  surface  to  which  we  apply  our  algorithm  to  perform  scattered  field 
and  power  computations.  In  what  follows  we  list  the  major  tasks 
accomplished  and  the  report  in  which  details  were  given. 

(a)  Generation  of  a  statistically  known  random  surface 
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This  task  was  reported  in  the  quarterly  report  for  the  period,  August  to 
October  1987.  It  was  shown  that  one  can  specify  the  surface  height 
statistics  and  its  correlation  function  and  that  the  algorithm  would  then 
generate  surface  patches  with  the  specified  statistical  properties. 

(b)  Specification  of  minimum  simulation  parameters 
The  minimum  number  of  surface  points  per  electrical  wavelength  or  per 
correlation  length  (whichever  is  the  smaller)  was  reported  to  be  five  in  the 
quarterly  report  for  the  period  August  to  October  1987.  The  minimum 
effective  width  of  the  antenna  pattern  was  found  to  be  at  least  three  times  the 
surface  correlation  length.  One  must  integrate  the  antenna  pattern  until  it 
dropped  to  10*3  of  its  value  along  the  boresight.  Additional  studies  on 
simulation  parameters  were  reported  in  the  quarterly  report  for  the  period 
November  to  January  1988.  There  it  was  indicated  that  to  obtain  an 
acceptable  estimate  of  the  scattering  coefficient  one  needed  to  average  at  least 
25  samples.  It  was  also  found  that  the  minimum  width  of  the  illuminated 
area  should  be  either  seven  surface  correlation  lengths  or  seven  electrical 
wavelengths— whichever  was  the  larger. 

(c)  Integral  equation  method 

The  issue  of  rough  surface  scattering  depends  on  obtaining  an  accurate 
representation  of  the  surface  current.  Once  the  surface  current  is  known 
only  a  two-dimensional  integration  has  to  be  performed  to  arrive  at  the 
scattered  field.  It  is  well  known  that  the  surface  current  is  defined  by  an 
integral  equation.  Hence,  an  attempt  was  made  to  solve  this  equation  and 
arrive  at  an  estimate  for  the  surface  current  (see  quarterly  report  for  the 
period  February  to  May  1989).  We  shall  refer  to  this  approach  as  the 
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integral  equation  method.  To  verify  the  applicability  of  this  estimate  it  is 
best  to  compare  scattering  coefficients  based  on  it  with  those  based  on  the 
moment  method.  Since  in  practice  we  cannot  use  the  moment  method  in  a 
three-dimensional  surface  scattering  problem  due  to  its  excessive  matrix 
size,  we  have  performed  scattering  calculations  for  two-dimensional 
scattering  problems.  This  was  done  and  reported  in  the  quarterly  report  for 
the  period  February  May  1989.  The  comparisons  between  the  two 
approaches  indicate  that  except  for  surfaces  with  rms  slopes  larger  than  0.4 
the  scattering  coefficient  obtained  using  the  integral  equation  method  is 
accurate  at  all  frequencies  and  for  all  incidence  angles  from  nadir  to  60 
degrees  from  the  vertical.  In  view  of  this  success  the  same  concept  to  obtain 
a  current  estimate  is  used  in  three-dimensional  surface  scattering  problems. 
It  was  shown  in  our  quarterly  report  for  the  period  February  to  May  1989 
that  simulated  results  could  be  obtained  and  compared  with  measurements  as 
well  as  theoretical  predictions.  The  only  drawback  with  the  integral  equation 
method  is  that  it  takes  considerable  amount  of  computer  time  to  do  the  three- 
dimensional  surface  scattering  problem  because  one  has  to  evaluate  a  four¬ 
fold  integral. 


(d)  Frequency  domain  method 

In  an  effort  to  find  a  more  efficient  approach  we  have  developed  a  frequency 
domain  method  (described  in  the  quarterly  report  for  the  period  from  May  to 
August  1989)  where  we  use  the  standard  spectral  representation  for  the 
scattered  field  and  the  boundary  condition  as  the  integral  equation  to  be 
solved.  When  applied  to  a  two-dimensional  surface  scattering  problem 
considerable  saving  was  obtained:  what  takes  an  hour  on  a  VAX  8700 
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computer  using  moment  method  was  reduced  to  3  minutes.  Efforts  are 
currently  underway  to  extend  this  method  to  three-dimensional  problems. 

3.  FREQUENCY  DOMAIN  METHOD  IN  3-D  SCATTERING 

The  theoretical  basis  for  the  frequency  domain  technique  in  a  three- 
dimensional  scattering  problem  is  given  below.  The  concept  is  the  same  as 
in  a  two-dimensional  problem.  Thus,  the  form  of  the  scattered  field  is  taken 
to  be 


Es  = 


4 

dky  a(kx,  ky)  e-ik‘x+Jkyy  -Jk*  • z 


and  it  should  satisfy  the  tangential  boundary  condition,  nx  E  =  0.  For  a  y- 
polarized  field  incident  in  the  xz-plane,  the  x-component  of  the  tangential 
field  is 


Ej,  +  Ey  +  Zy  E!  =  0 

Since  surface  depolarization  is  usually  more  than  an  order  of  magnitude 
smaller  than  the  like  and  surface  slope  is  usually  less  than  unity,  for  like 
polarization  we  can  drop  the  third  term  in  the  above  equation  and  write 


g-jkwX  +jkjoZ(x,  y}_  _ 


a  (kx,  ky)  e-i^+j^y  -ik*  • z  (*•  y)  dkxdky 


(1) 
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In  (1)  the  expression  on  the  left  hand  side  is  the  incident  field  and  a  is  the 
unknown  to  be  determined.  To  do  so  we  multiply  both  sides  of  (l)by 

e  jkl°X  jkyoy  +  jk2°Z^x’  and  the  pattern  function,  G  (x,y)  and  integrate  over  the 
illuminated  area  yielding 


0-j2kIox-jkyoy  +  j2kzo  z  (x,  y)Q  yj  (jxdy 


ajkx,  k  y)G  (x>  y)  dxdy  dk^  dk 


(2) 


Let 


I(kx,ky)=  f  I  ejik--k”)x  +  jfk>-k-)y-jfk‘-k”)z,x-y,G(x,  y) dxdy 

Note  that  I  (kx,  ky)  is  a  function  which  has  a  sharp  peak  around  kox  and  k 
Hence,  we  can  rewrite  the  right  hand  side  of  (2)  as 


ov. 


! 


and 


(kx’ky)I(kx-ky)dkxdk' 


=  a(k  ,  k  ) 

l  xo’  yo/ 


I  (k  .  k  )dk  dk 

\  x'  y/  x  y 

(3) 


a 


g-J^xoH^+j^x.y) 


G(x,y)dxdy 


where  I  a 

O  “ 


Ifk  ,  k  j  dk  dk 

\  x  yl  X  y 


(4) 


Next,  we  need  to  find  the  relationship  between  the  scattering  coefficient  and 

a. 
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To  do  so,  let  us  consider  the  ensemble  average  of  the  following  product 


a(kx’  ky)a(kV  k'y)  )  =  (f 


^-jzkjx-x')  +jz(kx-k't)  x'  -  jzky(y  -  y')  +  jz  (k,  -  k'y) 


e-4o2{kl2  +  k'z2-2kzk'z  p)/  2q  |x?  y|  q  (x',y')  dxdy'  dydy' 


-jzkx  %  +  jz  A  kx  x'  -  jzky  ti  +  jz  Aky  y' 


3  2<T^k  +  kz  2kzkl  P  ^G  (£+x',  T)  +  y')  G  (x',y')  dx'dy'  d^dr|' 


where  p  is  the  surface  correlation  function  and  a  is  the  surface  rms  height. 
If  k  =  k'  ,  k  =  k'  ,  k  =  k'  ,  we  obtain  another  expression, 

x  x  y  y  z  z  r 


a(k,'ky)a’(k,"ki  =  {|a(kx'ky[ 


;-jzk,(x-x')-i2k,(y-J,-)-(2k,<,)!(l.p)G(!Ci  y)G(x.  y-)dxdxdyd/ 
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<1  ■ -Mg  (C  +  x-,  „  +  y'j G (x',  /)  dyndx'  dy 

(6) 

Let  us  assume  that  the  pattern  function  is  Gaussian, 

G  (x,  y)  =  e-^Vi/g2 

Then,  the  following  integrals  given  below  by  (7)  and  (8)  can  be  evaluated 

^  ,  f  f  (x+C)2  x2  (y  +  T|)2  y2 

g  g  g"  g‘dxdv 


G  (x+^,  y+T))  G  (x,  y)  dxdy  = 


=  e  2g  I  e 


-2K)/8dxe-jf  e"2(y  +  2)/8 


dy 


=  c  2g 


(8Vf )=Te  282 


(7) 


ejzAkiX+  jZ^k^G  |^  +  X?  T|+yj  Q  ^  yj  ^^y 


~  g2^  p 

l  £ 

=  - —  e 

2g  2  . 

e 

[2s2  •  2  j 

(8) 


Using  (7)  and  (8)  we  can  rewrite  (5)  and  (6)  as  (9)  and  (10). 

a-jzk  ^  -  j2k  rj  -  2  ajk*  +  k\  -  2kk\p(q,  ti  )) 
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(9) 


2  2  2  2  2  2 
e  2g2  ‘  2  2«2  y  2  d^dTi 


2g 


2  2 
C,  +T) 


e-j2k,?-j2k(,-(2k,0)(l-ri4.1))e 

o 


(10) 


Now  we  are  ready  to  calculate  the  ensemble  average  of  the  scattered  field 
with  the  aim  of  expressing  it  in  terms  of  (10)  as  follows 


Ef)=j  J  J  J  (a  (kx>  ky)  oc  (k'x,  k'yf)  ei(k*-k'*)  x+)(ky-k'y)-  j(k  -k'*)zdkxdk'xdkydk\ 

(11) 

Substitute  (9)  into  (11)  and  integrate  kx'  and  ky'  out  ,we  have 


pi,wV«. 


2  g 

a2 


q)2  ki^  -  j2kyT|-(2kz  of  ( 1— p(C.  T|)) 


exp  ( -  (42+  r|2)/  2  g2]  dij  dr)  dk,  dky 
Substituting  (10)  into  (12)  we  have 


(12) 


0\  f 

)  _2k 

E 

'  2 

s  J 

°iK 


(13) 


From  Ulaby  et  al  (1982)  Chapter  12,  the  scattering  coefficient  is  related  to 
the  integrand  of  (13)  as  follows 


a°  =  47tks2cos2es^(ja|kx,  kyf 


(14) 
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Since  a  is  given  by  (4)  we  can  evaluate  the  scattering  coefficient. 


4.  CONCLUSIONS 

We  have  accomplished  all  the  tasks  proposed  and  gone  beyond  them  by 
investigating  new  approaches  as  shown  in  the  latter  part  of  Section  2  and 
Section  3. 
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